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Abstract 

We present a simple proof of monotonicity for cooperative systems of 
linear ordinary differential equations, without having recourse to approx- 



\^ I imation by strongly cooperative systems. 

i-G . A system of linear ordinary differential equations (ODEs) 

x' = A{t)x, t€{a,b), (1) 



where A{t) — [aij{t)]f,-^ is a matrix with real entries, is called cooperative if 
for each t £ (a, b) the off-diagonal entries of the matrix A{t) are nonnegative. 
^ ' We assume in the remainder that the matrix function A: {a,b) — > R"^" is 

^J^ , continuous, so for any {tQ,XQ) G (a, 6) x M" there exists a unique solution 

\q'. [{a,b) Btt^ x{t;to,xo) e 



of the initial value problem 

jx' = A{t)x, 






/^ . An important property of cooperative systems of ODEs is that their solution 



operators have some monotonicity properties. To formulate them, write 

R'J: :^ { X = (a;i, . . . , x„) g M" : x, > for all 1 < « < n }, 
M"+ := { X = (a;i, . . . , x„) G M" : Xj > for all 1 < « < n }. 

(Ml) For any <o £ {a,b), if Xq G K" then x{t;to,xo) G M" for all t G {to,b). 
(M2) For any to G (a, 6), i/ Xq G K"_|_ then x{t; to, Xq) G K"_|_ for all t G (to, ^)- 

The approach to proving the monotonicity (more precisely, property (Ml)) 
in all the literature on cooperative systems I have consulted (cf., e.g., [2], [4], 
[3], [5]) rests on the following. 



• A cooperative system is approximated by systems whose matrices have 
their off-diagonal entries positive (strongly cooperative systems). 

• It is proved that in a strongly cooperative system, if a nonzero xq is in 
K" then x{t;to,xo) G ^++ for t G {h,b)- That property is called strong 
monotonicity. 

• With the help of the continuous dependence of a solution on a parameter 
one obtains (Ml) for the original system. 

Indeed, the above approach is usually applied to general nonlinear systems, and 
from that point of view linear systems are special cases of nonlinear ones. 

When one wants to obtain (the apparently stronger) property (M2) one 
needs to apply some additional ideas, for instance, that the solution mapping is 
a homeomorphism (for a different method, see [2])- 

Below I present a simple and elementary approach, to the best of my knowl- 
edge (strangely) overlooked so far. 

Proof of (M2). Fix <o e (a, b), xq e M'J:_|_, and write x{t) = {xi{t),..., a;„(t)) 
for x{t;tQ,xii). We have 

i i j^i i ^ k j^i 

> ^(a^,(t)a;^^^J) =trA(t)]^x,. 

i j^i i 

So, denoting £^{t) :— xi{t) ■ . . . ■ x„(t), there holds 

e{t)>tTA{t)m, te{a,b), 

as well as £,{to) > 0. By continuity, £^{t) > for t > to sufficiently close to 
to- Suppose to the contrary that there are 6 G {to,b) at which £^{9) = 0. Let 
T G (to, b) be the least such 9. Since, for each t G (to, r) there holds 



C(0>^(^o)exp(/'trA(s)ds) 



to 

by letting t /^ t we obtain a contradiction, as the left-hand side tends to zero 
whereas the right-hand side goes to a positive number. 

It suffices now to notice that this means that Xi{t) > for all t G (to, b) and 
all indices iG{l,...,n}. D 

It should be mentioned that I came across the idea of using the product of 
the coordinates while reading the paper p"' by M. Benai'm and S. J. Schreiber. 



1. As the mathematical tools used are rather elementary, the proof is acces- 
sible to students who have taken a course in elementary calculus, linear 
algebra and differential equations. 

2. The proof carries over, almost word by word, to the case of Caratheodory 
conditions. 

3. It would be interesting to give an analog of it for the case of systems of 
linear ODEs preserving an arbitrary cone in R", compare e.g. [S]. 
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